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Abstract 

We analyse a generalised Gross-Pitaevskii equation involving a 
paraboloidal trap potential in D space dimensions and generalised 
to a nonlinearity of order 2n + 1. For attractive coupling constants 
collapse of the particle density occurs for Dn > 2 and typically to a 
5-function centered at the origin of the trap. By introducing a new dy- 
namical variable for the spherically symmetric solutions we show that 
all such solutions are self-similar close to the center of the trap. Exact 
self-similar solutions occur if, and only if, Dn = 2, and for this case of 
Dn = 2 we exhibit an exact but rather special D = 1 analytical self- 
similar solution collapsing to a 5-function which however recovers and 
collapses periodically, while the ordinary G-P equation in 2 space di- 
mensions also has a special solution with periodic ^-function collapses 
and revivals of the density. The relevance of these various results to 
attractive Bose-Einstein condensation in spherically symmetric traps 
is discussed. 
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The experimental discovery of Bose-Einstein condensation (BEC) in trapped 
vapours of cooled alkali atoms ||, ||, |] has opened up unique possibili- 
ties for the investigation of collective many-body effects in dilute gases. In 
the experiments the cloud of atoms is isolated from the environment by a 
magnetic trap. After cooling the cloud exhibits Bose-Einstein condensation 
i.e. the existence of a macroscopically populated quantum state. The study 
of the dynamics of this quantum state is an important fundamental prob- 
lem in many-body quantum physics. For three space dimensions D = 3 
the dynamics of the condensate can be described within the Hartree-Fock 
approximation by the Gross-Pitaevskii equation 

+ i^A,<l> - - = 0, (1) 

2m m 

where ^(x, t) is the wave function of the condensate, the external potential 
V{x) models the wall-less confinement (the trap), m is the mass of an indi- 
vidual atom, is the scattering length, and A^, = J^t is the Laplace oper- 

i 

ator. A convenient choice for the confining trap is the paraboloidal potential 

2 

assumed here to be spherically symmetric for simplicity, i.e. V = ^^^x^. 
In this paper we are concerned with condensates in D = 3 and D = 2 di- 
mensions. The Bose-Einstein condensate in two space dimensions is only 
marginally stable in that below the critical temperature correlations decay, 
but decay only as a power law f^, ^. Recent experimental techniques allow 
realisation of a two-dimensional trap for e.g. spin-polarized hydrogen ad- 
sorbed on a helium surface [|^, If. The dynamics of trapped Bose-Einstein 
condensates and the search for the related soliton-like solutions of the Gross- 
Pitaevskii equations is thus an interesting and relevant problem also in two 
dimensions. In this paper we concentrate on some aspects of this dynamics 
and on the existence of self-similar solutions of the Gross-Pitaevskii equation 
in particular. Self-similarity is an important and useful concept in nonlin- 
ear dynamics, particularly so when collapsing systems are being considered 
|10[ as they are below. This phenomenon of collapse appears in Bose- 
Einstein condensates with negative scattering length, as for example in ^Li 
(see e.g. fill). In this paper we show that self-similar behaviour only appears 
in two-dimensional traps although 'attractive' condensates (as < 0) collapse 
for all D>2. 

To begin with we consider a generalized D-dimensional Gross-Pitaevskii 
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equation, which for units such that h = 1, m = 1/2 can be expressed in 
the form 

+ A^?/'-2kV|V'P"-^^¥ = 0. (2) 

Here A^,. is the D-dimensional Laplace operator and = Y.F xf. Notice that 
" generahsation" means here an exponent 2n instead of the 2 which appears 
in the ordinary G-P equation. 

We consider only the attractive case of Eq. (|^) k < and the boundary 
conditions are vanishing at infinity. An observation is that a symmetry which 
leaves Eq. invariant is 



in which rf^ is an arbitrary vector in D dimensions, (fo is an arbitrary phase. 
This symmetry reveals the, in general, oscillatory character of the wave 
packet dynamics of Eq. whether k > or k < 0. In Ref. and its 
references 'collapse' was demonstrated for uo = Q and k < 0. Solutions be- 
come singular in a final time interval if the condition 

nD>2 (4) 

is fulfilled. We show here how the same condition arises in the present con- 



text, where 7^ (and k < 0). Ref. has addressed the same problem 



(of 7^ 0) for n = 1 and D = 2 and D = 3. Following both and we 
use the functional U[iIj] = Jj^d r^\'ip\'^d^x in which r = : U[ip] > 0. From 
Eq. (^) this functional satisfies a second order ordinary differential equation 
whose solution is 

,1 4sin^(co't) ^ 9/ X T sin(2a;t) 



^ ' 



t 



with 



t/0 = t/M|t=0, Jo = ^f/M|*=0, ^iVL5 = - ^^0, (6) 
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where 



Jr° \ n + 1 4 / 



is an obvious "energy" functional and is the Hamihonian of Eq. (|^) with the 
bracket {ip{x), ip*{y)} = i6{x — y). 

Hamiltonian Eq. (0) is a constant of the motion fixed by the initial data. For 
K < and smooth enough initial data it is not bounded below while E^ls 
as defined in Eqs. (P) has the same properties. The condition E^ls < 0, 
for example, still admits a large amount of physically accessible initial data. 
A second constant of the motion is /^d \ip\'^d^x = A/", the total number of 
bosons (atoms). Careful scrutiny of U[ip] of Eq. then shows (see also [^], 
[|1^) that provided that 

K<0, Dn>2, Enls < 0, (8) 

with the exception of the special case Dn = 2, Enls = = 0, there is 
always at least one point t = G (o, ^ such that the right hand side of 
Eq. (^) becomes negative for t > t^. Since by its definition the functional 
U[^p\ is nonnegative, this contradiction leads to the conclusion that ip cannot 
be continued beyond the point t = and must exhibit a singularity. We 
show below that this singularity is typically N'5{x). However, for 

the special case Dn = 2, Enls = = 0, the functional U[ip] = f/o cos^(a;t) 
never becomes negative. We show below that collapse in l^/^p occurs with 
IV'P — >■ MSi^x) as t — > t*, but now this can be followed by revival and periodic 
collapse of period t^/uj. There is some evidence that a form of collapse could 
occur in general even when U[^] apparently remains positive, i.e. at some 
point t < t^: (see fl^, ^ and references therein where uj = 0). We shall 
assume here that collapse occurs only at a zero of U[iIj]. 
Thus the conditions Eq. (§) are sufficient for U[ip] to reach a zero at t = 
^* < ^ and, generically at least, {ipl"^ — >■ MS{x) there. These conditions 
are sufficient but not necessary: for given such evolution for Enls < 0, the 
transformation Eq. (^ can increase Enls to > while the evolution remains 
singular. This is true for example for the exact analytical solution Eq. (p5|) 
for Dn = 2 we give below. The formation of these singularities may be very 
sensitive to the initial conditions and the values of the parameters. Evidently 
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these results mean that for E^ls < initially collapse and blow-up will occur 
for all M > Mc (see also |13| and references therein). 



We turn to the problem of similarity solutions which within the terms of our 
analysis arise only for Dn = 2. We seek spherically symmetric solutions of 
Eq. (H) in the form '?/'(r, t) = A(r,t)e^'^'^^''^^ in which r = \x\. From Eq. (j^) we 
arrive at the set of equations 

= (9) 




It is not evident how similarity solutions could be constructed from this set 
of equations in the general case, and we therefore choose to make an ansatz 
for the amplitude variable A{r,t): 

Mr,)^(^f{;^fAM (U) 

This ansatz solves Eq. provided that the function rj satisfies 

| + 2f^^ = 0. (12) 

at or or 



Notice that the function Ao(?7) is arbitrary and the ansatz Eq.(|TT]) describes 
an arbitrary spherically symmetric solution. The gradient |^ is related to 
the velocity of the particles of the condensate, and, through Eq. (|12[), ri{r,t) 
is then related to the local time dependent concentration of condensate par- 
ticles. In fact ri{r,t) completely determines this concentration as is evident 
from the number of particles n{r,t) in the interval [0,r], which is 

n(r,t) = Qn / r''''A\r,t)dr = Qn / v''-'Al{r])dr^, Qn = j^rj^y (13) 
U^/^J 

while n{oo,t) = A/" is independent of t. From the ansatz Eq. (pH]) we can 
deduce that ri{r,t) is a monotonically increasing function of r i.e. |^ > 0, 
and 7] —>■ oo when r oo. Also, in the vicinity of the origin r = 0, r] behaves 
as 

ri{r,t) =r/p{t) + 0{r^), (14) 
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where p{t) is a function of time. The solution for r]{r,t) is self-similar if 
7] = r/p{t) exactly. This 'self-similarity' is in the sense that the function 
f] depends now on a single variable t] = r/p{t). From Eq. ( ]T^ ) it follows 
immediately that in this case the phase (f){r, t) is quadratic in r, 

0(r,t) = 0o(t) + ^^rl (15) 

The Eq. (|10D should now be understood as an equation for Aq. Consider first 
the case nD = 2. Separating the variables in this equation we find that 

^1 (^""^) " ~ + = ' ^''^ 

<Po + y,=0 (17) 



p" + ^V-^ = 0. (18) 



4A 

Here A and p are arbitrary constants. 

A solution of Eq. ( |18|) can easily be found in the form 



p(t) = ^/cos2(u;t) + ^ sin^(u;t). (19) 



Other solutions can be obtained through the transformation t ^ t + to and 
p(t) — * h(t)p{s(t)), where 



h{t) = {Vl + + acos{2ujt))^, s{t) = -tan"^ ((yiTa^ - «) tan( 

(20) 

We have thus demonstrated that for nD = 2, r]{r,t) = r/p{t) with p{t) 
given by Eqs. (|19D and (|20|) is indeed a solution, and Eqs. ([TT|) and (|l^) 
now provide the corresponding self-similar solution of the Gross-Pitaevskii 
equation Eq. (0). For the explicit form of this solution one still needs to solve 
Eq. (|16D for AqIv)- 

In the case nD ^ 2 there are no self-similar solutions (except the trivial case 
p =const). Indeed, for the existence of such solutions we need to require that 
both Aq" and Af^Ao/AQ are functions quadratic in rj. These conditions cannot 
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obviously be satisfied. This means that even though the solution given by 
Eq. ( |TlD is locally self-similar for any D in the vicinity of r = 0, the exact 
self-similarity is only realised for Dn = 2. 

For the self-similar solutions there are two integral identities. Multiplying 
Eq. ([I6|) by ti^'^Aq and by rj^dAQ/dr], respectively, and integrating by parts, 
we find after a little algebra that 

Using the identity Eq. (|2T|) we easily find that the total energy of the solution, 
E[ip], Eq. (|^) is given by 



Em = -Ap) I r]''^'Al{r^)dr^, e(p) = (p')' + + ^- (23) 



As an example of an exact solution of Eq. (^) we consider here the attractive 
generalised G-P equation in one dimension: D = l,n = 2, X = and k < 0. 
In this case we find that 

MV)= I /° , P=l\'^\pi (24) 

'cosh [^JQ\k.\Pov) 

and the solution of Eq. (|]) can be expressed in the form 

^(^, t) = ^xp{-ztanM)(fx^-|;|^bg)} ^^^^ 

TJt) 



^os{u;t) ^cosh(|v^pg^ 



For an attractive condensate (k < 0) we expect the solution to become 
singular at a finite time. But it is indeed obvious that the solution Eq. (B^) 



becomes singular for t ^ when its amplitude diverges as 1/y ^ — t. In 
this limit I'j/'p |y^2^5(x) = Af6{x) which is the convergence to the S- 
function expected. Notice that j-j/'p from Eq. (B^) is now periodic of period 
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- while the solution Eq. ( pSj) itself has the jumps in phase, compounded 
by branch point singularities, when crossing the singularities of \ip\ at t = 
^{2k + l), k G Z|§. 



If now we simply assume that the point of collapse is t = defined below 
Eq. (I), it can still be shown that the collapse occurs to a (5-function centered 
on the trap. A consideration leading to this conclusion is the following: the 
equality U[ip] = means that \ip{x,t^)\ = for any x except possibly at the 
origin. Since Af = J d^xlipl"^ is a constant of motion identified as the total 
number of bosons (atoms), the obvious physical solution is {ipl"^ = AfS{x) 
excluding other possible generalised functions. The spherically symmetric 
case can be treated rigorously. Consider for this the functional U[ip] taken 
on the ansatz Eq.(|rTp, i.e. 

oo 

U[^]= j dr,r\r,,t)r,''-^Al{r,). (26) 



For U[^] = it immediately follows from Eq. (|26D that r{ri,t^) = 0. For an 
appropriate arbitrary test function (p{r) consider now the limit 

/•oo 

\imnn / drr^~^\ip(r,t)\'^(p(r) 
t-^t, Jo 

= \imnn dr]r]''-'\Ao{v)\Mr{v,t))=^viO). (27) 
*^** Jo 

This result means rigorously that for spherically symmetric solutions for 
which U[ip] evolves to a zero at t = t*. the system 'blows-up' to the S- 
function singularity 

lim \^{r,t)\^ = J\f5{r). 

This result is particularly evident for the self-similar solutions for which rj = 
r/p. In this case 

/ A \ \ /•rvn 

D+l a2i 



and the point of collapse for A < {E^ls < 0) can be readily found as 

U = (l/cu)tan-i (uj/2^\). 

Notice again that when A = the functional U never becomes negative 
and there is a possibility of periodic (5-function collapses and revivals of the 
condensate density in this case of two-dimensional traps. 
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Even though the methods are different, some part of the results reported here 



is analogous to that obtained in Refs. |1J, [T^ for the Nonlinear Schrodinger 
equation (NLS), which is the Gross-Pitaevskii equation with u = 0. This 
analogy is related to the fact that for Dn = 2 the generalised NLS and GP 



equations are equivalent. For the change of variables ||T7 



e = - tan(cjt), Zi = J . (28) 

UJ COS{Ljt) 



D 



?/'(x,t) = (cos(cjt)) 2 expi^-i—tan{ujt)r^jp{z,6) (29) 
maps Eq.(^ to 

2k, 

ipe + A,p = 0, (30) 

and it is clear that for Dn = 2 the 6 dependence of the effective 'coupling 
constant' disappears and the NLS system is recovered. This means in partic- 
ular that the whole variety of results available for the two dimensional NLS 
for n = 1 is directly applicable to the Gross-Pitaevskii equation for n = 1 
in two space dimensions. It is interesting that the Gross-Pitaevskii equation 
only allows a self-similar solution of the type considered in this paper in this 
case when it can be exactly transformed to the NLS equation. 
It is worth mentioning that for Dn = 2 all self-similar solutions of the Gross- 
Pitaevskii equation are invariant under the transformation 

^.(M)^Mr^exp(||.^).(^..w). (31) 

For the solution Eq. (^) this transformation means a mere rescaling Pq — >■ 
Po/ia + Vl + a^)4. 

We emphasize that our similarity analysis of the Gross-Pitaevskii equation 
is based on the ansatz Eq. (p!T|). This approach is applicable to the Gross- 
Pitaevskii equation in D space dimensions and with an arbitrary external 
potential V{x). It can also be shown that the dynamics described by the 
Gross-Pitaevskii equation for an arbitrary initial condition which has an ex- 
tremum, is effectively equivalent to a system describing a Z)-dimensional 
classical particle. This dynamical system generalises that found in for 
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Gaussian initial profiles through a variational approach. These results will 
be reported in a forthcoming publication. 

We showed in this paper that for Dn = 2 alone the generalised Gross- 
Pitaevskii Eq. (0) equation allows self-similar solutions, and that in this case 
it can be exactly transformed to the NLS equation with no trap potential. 
An explicit solution was given for D = l, n = 2, k<0 which displayed a 
delta function divergence at a finite time. We further showed, for Dn > 2 
and K < 0, that all spherically symmetric solutions with Enls < collapse 
in a finite time to a (5-function centered at the origin of the trap while we 
showed generally that even without such symmetry evolution may be to the 
(5-function singularity. The ordinary Gross-Pitaevskii equation in 2 space di- 
mensions and with k < 0, Emls = 0, was shown to have periodic 5-function 
collapses and subsequent revivals of the particle density. 
One of us (AVR) wishes to thank M. Wadati, J. Hietarinta and S. Jaakkola 
with collaborators for useful discussions. 
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